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Since the first polar code of Arikan, the research field of polar codes has 
been continuously active. Improving the performance of finite-code-length 
polar codes is the central point of this field. In this paper, the parallel 
concatenated systematic turbo polar code (PCSTPC) model has been 
proposed to improve the polar codes performance in a finite-length regime. 
On the encoder side, two systematic polar encoders are used as constituent 
encoders. While on the decoder side, two single iteration soft-cancelation 
(SCAN) decoders are used as soft-in-soft-out (SISO) algorithms inside the 
iterative decoding algorithm of the parallel concatenated systematic turbo 
polar code (PCSTPC). As compared to the optimized turbo polar code with 
SCAN and BP decoders, the proposed model has about 0.2 dB and 0.48 dB 
gains at BER=10™), respectively, in addition to 0.1 dB, 0.31 dB, and 0.72 
dB gains over the TPC-SSCL32, TPC-SSCL16, and TPC-SSCL8 models, 
respectively. Moreover, the proposed model offers less complexity in 


comparison with other models, therefore requiring less memory and time 
resources. 


This is an open access article under the CC BY-SA license. 


Corresponding Author: 


Wallaa Yaseen Alebady 

Department of Electrical Engineering, College of Engineering, University of Babylon 
Babylon, Iraq 

Email: walaa.lafta@student.uobabylon.edu.iq 


1. INTRODUCTION 

Polar codes have changed the performance equation in the coding theory since the first polar-code 
of Arikan in 2009 [1]. The attention of most researchers in the coding theory field has been focused on the 
polar codes. The successive cancellation decoder (SC) is the first polar decoding algorithm that is proposed 
by Arikan in [1]. Since the polar code with the SC decoder is the first capacity-achieving code at the infinity 
code length. The performance of polar codes deteriorates at the finite code length and small SNR, this paper 
presents a turbo polar code to improve the performance of polar codes, especially at the finite-code-length 
and small SNR. There are several decoding algorithms suggested to improve the performance of polar codes 
at small and moderate code length. 

Tal and Vardy [2] propose the successive cancellation list (SCL) decoder. The SCL decoder offers 
excellent performance, especially when coupling the CRC code as a genie-code along with it [3]-[5]. At a 
long and moderate code length, the complexity of SCL is very high. Niu and Chen [6] propose the successive 
cancellation stack (SCS) to reduce the high complexity of the SCL decoder. The belief propagation (BP) 
decoder can also be used as a decoding algorithm for polar codes [7]. Fayyaz [8] proposes the soft successive 
cancellation (SCAN) decoder as a soft-in-soft-out (SISO) decoding algorithm. The other soft-in-soft-out 
algorithm is the soft successive cancellation list (SSCL) decoder, which offers a higher complexity as 
compared to the SCAN decoder [9]. Among the above decoding algorithms, only the BP, SCAN, and SSCL 
algorithms can be used as SISO algorithms with the turbo iterative decoding process. 
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Turbo-polar codes (TPCs) use the concatenation structure and the iterative decoding algorithm to 
improve the performance of the classic systematic polar codes (SPCs). The works of [9], [10] present turbo 
polar codes with SSCL and BP decoders as iterative decoding algorithms. Because of the large number of 
iterations and the list size of the BP and SSCL algorithms, respectively, the complexity of the previous turbo 
polar codes is very high. Therefore, the SCAN decoder is used as a SISO algorithm for the turbo iterative 
decoding algorithm [11], [12]. The results of these codes will be compared with the results of the present 
work later on in section 6. 

The remainder of this paper is organized as follows: section 2 describes the systematic polar 
encoding and the successive cancellation decoder. In section 3, the soft successive cancellation (SCAN) 
decoder is explained. The structure of the concatenated turbo polar code is presented in section 4. Section 5 
presents the complexity analysis of the presented turbo polar code and other codes. Finally, in sections 6 and 7, 
the simulation results and the conclusion are presented, respectively. 


2. PRELIMINARIES 
2.1. Systematic polar code encoder (SPC) 

To improve the bit error rate (BER) performance of the polar codes, Arikan proposes the SPC, 
which offers a BER performance over the non-systematic polar codes [13]-[15]. The SPC offers the ability to 
separate the parity bits and the information bits after the encoding process. The original polar coder can be 
used to create the systematic polar codes (SPCs) by modifying the original expression [13], [16]-[18], as (1): 


x=u*FO® =uxG 


F= h i = first order kronecker matrix (1) 


where G = F®"; u,x € FY, G € F% , and F®” is the nth Kronecker power of F. For designing polar codes, 
F must be fixed at the binary field. 

The binary vector u consists of two parts: u = (ure, Up) for ffc (1,2,...,N). The first part (upc) 
carries the user information, while the second part (uç) carries the frozen bits that are fixed in the encoding 
process and are known at the decoder. Thus, the encoded bits in (1) can be separated as: 


x = UpcGre + Ur Gf (2) 


where f = (1,2,3,...,N)\f°. Gre is a submatrix of G comprising of rows that correspond to indices f° and 
Gp are the same Gpe but with indices f. By selecting the length of the set f°, the code rate can be adjusted as 
R = K/N, where K is the length of set f°. For systematic polar codes (SPCs), the code in (2) can be 
rewritten as: 
X4 = Upc * Gee a + Uf Gf a (3) 
Xp = Ufc * Greg tur * Grp (4) 


where A and B correspond to indices of user and parity bits in SPC, respectively. Gre 4 is the sub-matrix of G 
comprising of elements (Gap) with a E€ f° and b E A, and likewise for the other sub-matrices. 

For Systematic Polar Codes (SPCs), x4 corresponds to upc and uf is fixed as frozen bits known as 
zero vectors. So, the use of systematic code can be calculated through (3) as (5): 


Ure = (x4 — ur Gr a) * Grt, = (xa) * G4 (5) 
by substituting (5) in (4), the parity bits (xg) of SPC can be calculated as (6). 

Xp = |(x4) * Grey] * Greg (6) 
From (5) and (6), the conclusion can be drawn that the systematic polar code can be achieved if (only if) 


f° and A have the same number of indices, and Gre 4 is an invertible matrix. Thus, (5) and (6) can be 
rewritten as (7) and (8). 
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upe = (xpe) * Grelpe (7) 
xs = [(xpe) * Gi ye * Gre (8) 


2.2. Successive cancelation (SC) decoder 

As the native decoding algorithm, Arikan proposed the SC decoder for the polar codes [1]. For polar 
code (N, K), a decoding tree of an n-depth can be used to illustrate the decoding process, where n = log,N. 
Figures 1 and 2 show the decoder structure and 3-depth tree for polar code (8,4), respectively. In decoding 
the tree and at t-stage, there are 2"* nodes. Each node receives and propagates afand ßt vectors, 
respectively, from and to the parent node, where 0 <t <n. The received soft information æf and the 
propagated hard information Bf consist of 2‘ values, as shown in Figure 2 (e.g. the first node of the second 
stage receives four LLRs a2, [a2,a?,a2,a2] from the root node and propagates four binary bits pĝ, 
[B3, B?, BZ, B2]). To compute the LLRs ætt of the left child of the t-node, the following expression must be 
applied: 


att = f (af, at, a) ) 


while the LLRs a{~* of the right child are computed after estimating the hard information B/~+ of the left 
child as (10): 


att = glat sat pa) = 


where l is a vector index ranging from to 2°. 
The message bits MÑ; are estimated based on the values of LLRs @?, where these LLRs can be 
computed as (11): 


o Peoia Pano) 
ai Loo ear GD 


where R is the received log-likelihood ratio (LLRs) vector of the transmitted codewords through a channel. 
Based on the following rule, the message bits 77; can be computed as (12). 


l 


z 0 
A = {° ifai 20 (12) 
1 otherwise 


The estimated value of the frozen bit is M;=0, regardless of the value of the a?. The feedback hard 
information Bf can be computed by a linear combination of some decoded hard bits that are previously 
estimated as: 


Pi“ = Br? ® eee (13) 

Biye- = Broe--1 (14) 
where B? = M. The f and g functions are defined as (15) and (16) [19]. 

f (a, b) = sign(a * b) + min(Jal, ||) (15) 

g(a, b, B) = [(-1)?*F xa] +b (16) 

For the systematic polar codes, the decoding process needs an extra step to obtain the systematic 
information bits. To estimate the information bits, the outcomes of the decoder (m) must be re-encoded (i.e., 
multiplied by generator matrix [G]) to obtain the systematic information ( BP). This extra step can be avoided 
by propagating the hard decisions (™,) in the decoding tree [20]. 

After the last iteration of the SC decoder, the resulting matrices œ and 6 contain information about 


noisy received LLRs, recovered encoded bits, and estimated message bits. B(:,1) contains hard reversed- 
indexed information for the recovered encoded bits, whereas a(:,n + 1) contains LLRs for message bits that 
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correspond to the hard decisions at B(:,n + 1). Noisy received LLRs at æ(:,1) are not updated throughout 
the decoding process. 
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Figure 2. Decoding tree of polar code with N=8 


3. SOFT-CANCELATION DECODER (SCAN) 

Based on the SC’s schedule, the SCAN decoder massively reduces the decoding complexity, as 
compared to the flooding schedule of the belief-propagation (BP) decoder [8]. In the SCAN algorithm, the 
soft information is propagated in both directions (i.e., from right to left and verse versa). Although the 
propagation of soft information in both directions somewhat increases the complexity and latency in 
comparison with the SC decoder, it improves the performance of the decoding process. The polar factor 
graph in Figure 1 of the SC decoder can be used with the SCAN decoder. In the beginning, the received noisy 
LLRs from the channel are put in vector œ? at left; where 0 < i < 2” — 1. The knowledge of frozen bits’ 
locations can be exploited by the SCAN decoder to initialize the vector pf by a priori information, where the 
frozen locations are filled by infinity because these locations correspond to zero bits. The locations of 
information may be zero or one equiprobably, so they are filled with zero LLRs as (17) [21]: 


pr ={” ifie f 


17 
0 otherwise ee 
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while computing qj with @j.;_, Which has already been computed; æf (where 0 < t < n,0 < j < 2” — 1) is 
partially populated from top left to bottom right. This partially computed @ matrix can be used as a priori 
information to update 6 from top right to bottom left. At the final step with fully populated a; (when t = 0 
and j = 2” — 1), all nodes of the factor graph contain soft extrinsic information (LLRs) which correspond to 
the values of the fully populated # matrix. At the first iteration, all values of p are zeros except př, 
where 0 <i < 2" — 1, so there is no extrinsic information that can be used by the SCAN decoder. At the 
second iteration and other successive iterations, the SCAN decoder exploits the extrinsic information 
computed with previous iterations [8], [22]. 

The updating process can be described as follows based on the factor graph of the basic polar kernel 
in Figure 3 [23]: 


aft = fab af + BED, af = af + f(a, 65) 


Bo = f(Bo *) Bi + ai), Br = Br’ + FBO", ao) 


where f is the box-plus operator (H) as (18) [24], [25]. 


xfy+ a = 2tanh"! [tanh (=) x tanh B) (18) 


The previous form of box-plus requires more memory resources and execution time, so a hardware-friendly 
approximation can be used instead of the previous equation as (19). 


f(x,y) =x By = sign * y) * min(|xl, lyl) (19) 
From Figure 4, it can be stated that the systematic polar codes with the SCAN decoder have the same frame- 


error rate (FER) as the non-systematic polar codes. In terms of bit-error-rate (BER), the systematic polar 
codes offer a significant advantage over non-systematic codes, as shown in Figure 4. 


10 T T T 
=E SPC-BER-(N=64,lo=3) 
annn E} SPC-FER-(N=64,10=3) 
ji t-1 a | i H 
ah Bi a i 0 ie i —P- NSPC-BER-(N=64,lo=3) 
etme “> NSPC-FER-(N=64,l0=3) 
> < > < ih 102 
a 
“ 
W 10° 
at Bi a B Wig 10% 
> + > + 10° 1 1 L 1 
25 3 3.5 4 45 5 
Eb/No (dB) 
Figure 3. The factor graph of basic polar kernel Figure 4. BER and FER comparisons between 


systematic and non-systematic PC (64,32) 


4. TURBO CONCATENATED POLAR CODES 

The construction of the turbo codes depends on several methods. However, the presented turbo polar 
code uses the traditional way that depends on two recursive and systematic codes at the encoder side. While 
at the decoder side, two iterative soft-in soft-out (SISO) decoders (such as SCAN, BP, and SSCL) are used. 


4.1. Turbo encoder structure 

Figure 5 shows the encoder’s construction of turbo polar code. This construction includes two SPC 
encoders: the first encoder receives the original information bits (d) and generates the corresponding 
encoded bits, whereas the second encoder receives the interleaved data bits (d7) and generates the 
corresponding encoded bits. The output of the construction shown in Figure 5 consists of three branches 
multiplexed with each other. The first branch carries the original user bits (d), the second branch conveys the 
parity bits of the first encoder (p+), while the third line carries the parity bits of the second encoder (p°). The 
multiplexed output of the concatenated encoder will be c = [d pt p?], which is modulated by the BPSK 
modulator and transmitted through the AWGN channel. 
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Figure 5. The structure of turbo polar encoder 


4.2. Iterative decoding algorithm 
The encoded information after being transmitted over the AWGN channel is decoded by several 
algorithms, successive cancellation (SC), successive cancellation list (SCL), successive cancellation stack 

(SCS), belief propagation (BP), soft successive cancellation list (SSCL), and soft cancellation (SCAN). 

However, most of these algorithms cannot be applicable with turbo polar codes (TPCs) because of their hard 

outputs (1 or 0). Only three of these algorithms contain the soft-in-soft-out (SISO) property and can therefore 

be used with turbo-polar codes, namely: belief propagation (BP), soft successive cancellation list (SSCL), 
soft cancellation (SCAN). 

As shown in Figure 6, the SISO soft cancellation (SCAN) algorithm has been used with the proposed 
turbo polar code (TPC). The mechanism of the Iterative decoding algorithm can be summarized as follows: 

a) The received LLRs from the channel are separated into three vectors, namely, R4, Rp1, and Rp2, which 
are LLRs of the systematic bits, the parity bits of the first encoder and the parity bits of the second 
encoder, respectively. 

b) Ra, Rpi, and the scaled prior information (94) are the inputs of the SCAN-1 decoder that has one output 
(b,). The internal procedure of SCAN-1 is shown in Figure 6 (The small red box on the top right). 
Based on the (f°) and (f) sets, R4 is merged with Rp1 and used to initialize the first column of æ matrix 
of SCAN-1 decoder, where the vector (Ra) corresponds to f © indecies while the vector (Rp1) 
corresponds to f indecies. 


c) The extrinsic LLR (€,) can be calculated as follows [9], [11]: £4 = b4 — Aat 


o2 


pi, where o is the 
noise variance. 

d) The extrinsic LLR (€,) is interleaved and passed to the SCAN-2 decoder as prior information (P,). 

e) Ri, Rp2, and the scaled prior information (3) are the inputs of the SCAN-2 decoder that has one output 
(b2). The internal procedure of SCAN-2 is the same as SCAN-1 (see Figure 6: The red box). 

f) The extrinsic LLR (ez) can be calculated as follows [9], [11]: £, = b, — ae - 9%. 


g) The extrinsic LLR (e2) is de-interleaved and passed to of the SCAN-1 decoder as prior information (@, ). 


h) After the last outer-iteration, the estimated bits (d) are obtained by applying the hard decision to the 
output of the SCAN-2 decoder (bz) and then de-interleaving it. 
The stopping mechanism of the iterative decoding process can be developed based on two criteria. 
The first one depends on the number of outer iterations: after the last iteration, the iterative decoding process 
is finished. The second criterion depends on a particular rule. If this rule is satisfied, the iterative decoding 
process is terminated. Immediately after each iteration, the estimated vector (£) is calculated to be used for 
estimating the vector (û), as [12]: 


A = 2 [Gy] t = 2 * Gy 

The frozen bits of the vector (u) are fixed with polar codes and known at the receiver, so these bits 
with the frozen bits of the vector (&) can be used to satisfy the second criterion of the stopping mechanism 
through the: 

S= ulf) + af) 


if the S factor equals zero, the iterative decoding process stops. Otherwise, the decoding process continues [12]. 
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Figure 6. Iterative decoding of the proposed model 


5. SCAN ITERATIVE DECODING COMPLEXITY ANALYSIS 

As stated previously, only the SCAN, BP and SSCL decoders can be used with the turbo iterative 
decoding process; so, in this section, the complexity of these decoders is described and compared. A large 
amount of complexity is consumed by the SISO-1 and SISO-2 decoders, while the other components 
consume less complexity compared to the complexity of SISO-1 and SISO-2 decoders. As a result, it is 
possible to compute the complexity analysis of iterative decoding without taking into account the complexity 
of these components. 

The complexity of a single SCAN decoder is O(I; * NlogN); where I,,is the number of inner 
SCAN iterations. Thus, the complexity of a turbo iterative decoding based on a SCAN decoder is twice the 
complexity of a single SCAN decoder if the outer iteration equals one (louter = 1). The general formula of 
the SCAN iterative decoding process is: Iouter * (2 * Is, * NlogN). 

Table 1 shows the total complexity formulas of the iterative decoding with SCAN, BP and SSCL 
decoders. Isand Ig, represent the inner iterations numbers of the SCAN and BP decoders, respectively, whereas 
List is the list size of the soft successive cancellation list (SSCL) decoder. The BP decoder requires sixty inner 
iterations (i.e. Iz, = 60) to approach the performance of the SCAN decoder when working with only two inner 
iterations (i.e. /;, = 2). Obviously, the inner iterations of the SCAN decoder are smaller than the inner iterations 
of the BP decoder. For the SSCL decoder, the list size (Lyist) is always 32 or 64 [2]. It can therefore be stated 
that the complexity of turbo iterative decoding with SCAN decoder is smaller than the others. 


Table 1. Complexity formulas of different SISO decoders 


Decoder Scheme Complexity Formula 
SCAN Tourer * (2 * 15, * NlogN). 
BP louter i (2 j Ip; $ NlogN) 
SSCL louter id (2 i Liist ia NlogN) 


6. RESULTS AND DISCUSSION 

In this section, the simulation results (Figure 7) of the proposed model are presented in detail and 
compared with the results of other models shown in [9], [11], [14]. All results in this section are carried out 
over the AWGN channel and with a single-iteration SCAN algorithm. The maximum number of simulation 
frames varies according to the SNRs, with the upper limit of the frame errors being 1000 errors. 

Figure 7(a) shows the performance of the proposed scheme in terms of bit error rate (BER) and 
frame error rate (FER). These results are determined with the following parameters: i) Iouter = 6, ti) Iscan = 
1, iii) SF = [0.63], and iv) SF? = [0.75]. 

From Figure 7(a), it can be seen that the PCSTPC with a (128.64) SPC offers a more reliable 
performance compared to standalone (128.64) SPC with a SCAN algorithm. At BER=10~*, the proposed 


Turbo polar code based on soft-cancelation algorithm (Wallaa Yaseen Alebady) 


528 o ISSN: 2502-4752 


PCSTPC (128.64) has a 0.77 dB gain over the standalone (128.64) SPC. These results prove the theory 
proposed in [12], which states that “at finite code lengths, the performance of polar codes can be improved by 
using the turbo polar codes (TPCs)” [12]. 

Figure 7(b) compares the BER performances of the proposed model and the optimized model in 
[14]. The optimized model uses two decoder algorithms: SCAN and BP decoders. At BER=1074, the 
proposed scheme (with eight outer-iterations) has approximately 0.2 dB and 0.48 dB gains over the 
optimized models of SCAN and BP in [14]; respectively. At six outer-iterations, the BER-performance of the 
proposed and the optimized SCAN models approaches with each other. However, the proposed model offers 
an approximate gain of 0.35 dB over the optimal BP model (which has three outer-iterations and 60 inner- 
iterations). 


—©-— [14] C1, SFs: all 1.0 - SCAN 
—*— [14] C1, SFs: all 1.0 - BP 
—@- Proposed TPC(128,64)40=6 
=$ proposed TPC(128,64)-lo=8 
=E [14] C1, SFs:Optimal - SCAN 


—<- [14] C1, SFs:Optimal -BP 
10° 
10° | 
f ' 
i i 
4 | 
10 | $ 
' 
frn & 104} 
a D f k 
- ! 
10 i 
| 
H 
| 10° H 
=E PCTPC(128,64) BER H 
410° ©)» PCTPC(128,64) FER ; 
—©— PcTPc(256,128) BER i 
©" PCTPC(256,128) FER H 
—— Sys.SPC-SCAN(128,64) BER 10% i 
407 ||“ Sys.SPC-SCAN(128,64) FER | f 
' 
i 
3 
4 
10? L i 1 i J 
25 3 35 4 45 5 2.5 3 35 4 45 5 
Eb/No (dB) Eb/No (dB) 
(a) (b) 


Figure 7. Simulation results of (a) Comparison between the proposed PCSTPC and SPC-SCAN and (b) 
Comparison between the proposed PCSTPC and the optimized model [14] 


Figure 8(a) illustrates a comparison between the proposed model and the TPC-SSCL model 
(presented in [9]) in terms of BER performance. At BER=10~*, the proposed TPC offers 0.1 dB, 0.31 dB, 
and 0.72 dB gains over the TPC-SSCL32, TPC-SSCL16, and TPC-SSCL8 models, respectively. The 
proposed model with six outer-iterations also offers a good performance as compared to the TPC-SSCL 
models. 

In terms of FER, the performance of the proposed model is compared with the performance of the 
punctured turbo polar code (PTPC) model in [11] with the same (128.64) SPC component code. The 
proposed scheme has approximately 0.72 dB gain over the PTPC model, as shown in Figure 8(b). The 
parameters of all previous models are presented in Table 2. 

Finally, the proposed model has a superior advantage over all previous models in terms of 
complexity. The low complexity of the SCAN decoder and the low inner iterations (just single iteration) 
make the proposed model require less hardware and memory resources to be implemented practically. 
Although the TPC-SSCL32 model has a high degree of complexity, the proposed algorithm requires only 
eight outer-iterations instead of six to exceed the BER performance of the TPC-SSCL32 model, as shown in 
Figure 8(a). The mathematical expressions of complexity are shown in Table 2 for each model. These 
expressions are very crucial to determine the memory and time requirements for implementing any model. 
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Figure 8. Simulation results of (a) Comparison between the proposed model and the TPC-SSCL model [9] 
and (b) FER comparison between the proposed PCSTPC and the PTPC model in [11] 


Table 2. Summary of the studied models’ parameters 


Features [9] Model [14] Model [11] Model Proposed Model 
Enco. Comp. SPC and SPC SPC and SPC SPC and SPC SPC and SPC 
SISO-1 SSCL BP or SCAN SCAN SCAN 
SISO-2 SSCL BP or SCAN SCAN SCAN 
Interleaver Random Random Random LTE interleaver 
Puncturing x x v x 
IBP SCANorlist lust = 8, 16, 32, 64 Igp = 60, Iscan =? Iscan = 4 Iscan = 1 
Tourer 6 3,6 18 6,8 

(0) [2]outer Ol[2Touter O[2louter Ol2Iouter 
* (list * Nlog N)] * Upp or scan * Nlog N)] * (Iscan * Nlog N)] * (Iscan * Nlog N)] 
Complexity From Medium (with (BP>Very High) High Low 
lise = 8) to Very (SCAN>Medium) 


High (with List = 64) 


Iscan: Inner iterations number of SCAN decoder 
list: List size of SSCL decoder 

Touter: Number of iterations between the SISO 
algorithms 


Notes: 
Rp: Rate of SPC = K/N 
Tgp: Inner iterations number of BP decoder 


7. CONCLUSION 

The turbo polar code based on SCAN algorithm is introduced in this paper to improve the 
performance of the original polar code at a finite code length regime. The simulation results show that the 
presented model offers an excellent performance at finite-code lengths with low complexity as a result of 
using the single iteration SCAN algorithm. As a SISO decoder, the single iteration SCAN algorithm works 
perfectly with the iterative decoding algorithm, as each SCAN decoder helps the other to improve its 
performance. Therefore, the overall performance of the proposed TPC exceeds that of a standalone SPC with 
the SCAN decoder. From the simulation results, the conclusion can be drawn that the turbo polar codes 
(TPCs) have overcome the error floor problem, which is one of the major disadvantages in the turbo codes. 
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